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PART  I  (Introductory) 

•  Linear  IVPs,  Eigenvalue  problems,  linear  PDEs 

•  Manifolds  (“stay  on  manifold”  principle) 

•  Classical  problems  (“curved  path”  principle) 

PART  II  (Recent  results  on  exp  ints) 

•  A  unified  approach  to  exponential  integrators 

•  Order  theory 

•  Bounds  for  dimensions  of  involved  function  spaces 
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1.1  Linear  I  VPs 

One  may  for  instance  write 

ii  =  A(t)u ,  A  :  M  — ►  RnXn 

In  literature,  usually  u  e  Mn. 

LGI:  Magnus  series  or  related  (Cayley  etc) 

When/Why  use  this  scheme. 

1 .  Highly  oscillatory  ODEs,  large  imaginary  eigenvalues. 

Iserles 

2.  PDEs,  A(t)  unbounded,  classical  example:  Linear 
Schrodinger  equation  (LSE)).  Blanes  &  Moan, 
Hochbruck  &  Lubich. 

Recently  also  Landau-Lifschitz  equation  Sun,  Qin,  Ma 
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1.1  Magnus  works  on  LSE! 


i —  =  H(t )  u ,  H(t)  unbounded,  selfadjoint 
dt 

dex.pu  is  not  invertible  for  2k7zi  g  <t(u),  k  e  Z\{0}. 

Truncated  series  is  still  unbounded  at  oo. 

H  &  L  find  error  bounds  of  the  form 


C  hp  tm  max  II Dp  1w(£)|| 
o <t<tm 


D  is  a  “differentiation  operator”  related  to  the  LSE . 
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Eigenvalue  problems 


Stability  of  travelling  wave  solutions  to  PDEs.  Boils  down 
to  eigenvalue  problem 

Y  =  A(t,  A)  Y 

where  A  is  a  parameter. 

Needs  to  be  solved  for  several  A. 

Magnus  integrators  used  with  success  by  Malham,  Oliver 
and  others. 

Early  work  by  Moan  on  such  problems. 
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1.2  Problems  on  (nonlinear)  manifolds 

A  large  part  of  the  applications  I  know  involves  the 
orthogonal  group  which  acts  transitively  on  either  of 

•  The  orthogonal  group  itself  (or  its  tangent  bundle). 

•  Stiefel  manifold,  (n  x  p  matrices  with  orthonormal 
columns) 

•  The  n  —  1-sphere.  (Stiefel  with  p  =  1) 
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1.2  Orthogonal  group  problems 


Most  used  examples  are  on  n  =  3  (3D  rotations):  Free 
rigid  body,  spinning  top,. . . 

Scheme.  Most  LG  Is  work.  RKMK,  Crouch-Grossmann,. . . 
combined  with  all  possible  “coordinates”  exp,  Cayley, 
CCSK  etc. 

My  evaluation 

•  Most  Lie  group  integrators  do  little  else  for  you  than 
maintaining  orthogonality. 

•  Poor  long-time  behaviour. 

•  Hard  to  get  reversible  /  symplectic  schemes. 

•  There  are  exceptions  (Lewis  and  Simo,  Zanna  et 
al.)  but  these  LGIs  seem  expensive. 
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1.2  Stiefel  manifolds 


Some  applications  which  involve  computation  on  Stiefel 
manifolds 

•  Computation  of  Lyapunov  exponents 

•  Multivariate  data  analysis  (optimisation,  gradient 
flows) 

•  Neural  networks,  Independent  Component  Analysis 

Demands.  Maintain  orthonormality.  Inexpensive  stepping, 
cost  0(np 2)  per  step. 

Schemes.  Most  LGIs  work.  3KMK, 

Crouch-Grossmann,. . .  combined  with  all  possible 
“coordinates”  exp,  Cayley,  CCSK  etc.  Most  of  them 
can  be  implemented  in  0(np2)  ops  per  step,  but 
special  care  must  be  taken. 
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1.2  Stiefel  manifolds 


My  evaluation 

•  Lie  group  integrators  meet  requirements  specified 
in  literature 

•  Long-time  behaviour  has  not  been  an  issue. 

•  Overall  judgement:  Lie  group  integrators  are 
competitive,  if  not  superior  to  classical  integrators. 

Sources 

•  Dieci,  Van  Vleck  [schemes,  but  also  general 
viewpoints,  Lyapunov  exponents] 

•  Trendafilov.  [Multivariate  data  analysis] 

•  Celledoni,  Fiori. [Neural  nets,  ICA] 

•  LGIs  for  Stiefel,  Krogstad,  Celledoni  +  O 


Applications  of  Lie  group  integrators  and  exponential  schemes  -  p.  9/60 


1.2  Other  manifold  applications 


•  Certain  PDEs  whose  solution  evolve  on  (copies  of) 
S'”-1,  Lie  group  integrators  have  been  used. 

•  Some  special  types  of  manifolds,  e.g  described  by 
quadratic  invariants  like  oblique  manifold,  DelBuono, 
Lopez. 
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1.2  Conclusions 


•  There  are  several  manifold  applications  where  Lie 
group  integrators  now  represent  an  alternative  choice. 
Recent  research  have  caused  implementations  to  be 
much  less  expensive. 

•  Their  best  feature  is  that  they  preserve  the  manifold.  I 
have  seen  little  evidence  to  suggest  that  Lie  group 
integrators  is  natural  for  maintaining  additional 
geometric  structure.  Is  there  hope  for  improvement  on 
this  point? 

•  The  development  of  Lie  group  integrators  has  added 
important  insight  in  the  integration  of  DEs  on 
manifolds.  Understanding  of  numerics  has  become 
less  dependent  on  specific  coordinates  and 
embeddings. 
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1.2  Conclusions  (2) 


The  activity  on  Lie  group  integrators  has  caused  progress 
in  fields  related  to  geometric  integration: 

•  Computing  the  matrix  exponential 

•  Computing  highly  oscillatory  integrals 

•  Analysis  of  split-step  schemes 

•  Exponential  integrators 

•  Algebraic  structure  on  trees,  Hopf  algebras 

•  Computation  with  the  BCH-formula 
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1.3  Curved  path  principle 


•  Classical  numerical  ODE-solvers  progress  solution 
along  straight  lines. 

•  Lie  group  integrators  map  a  straight  line  in  some  other 
space  (Lie  algebra)  to  phase  space  through  a 
nonlinear  map. 

•  Allows  for  much  more  general  “movements”. 

Two  excellent  examples  provided  by  Munthe-Kaas. 

1 .  The  northern  light  equations. 

2.  PDEs  with  perturbation  terms  by  affine  action 
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1.3  Northern  light 
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1.3,  II  Time  Integrators  for  Nonlinear  PDEs 

Many  PDEs  are  of  the  (abstract)  form 

ut  =  L  u  +  N  ( u ) 

L\  unbounded  linear  operator  (like  A) 

N(u)\  a  (relatively)  small  nonlinear  term. 

Includes:  NLS,  Nonlinear  heat  equations,  KdV, 

Allen-Cahn,  Kuramoto-Sivashinsky,  and  many  more. 

Unbounded  L  requires  a  form  of  implicit  integrator. 

One  wants  an  explicit  scheme  for  the  nonlinear  part. 

Many  time  integrators  are  known  for  this  purpose. 
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Main  Classes  of  Time  Integrators 


LI  Linearly  implicit  or  IMEX  methods  (Varah,  Ruuth, 
Ascher,  Russo  and  many  more) 

SS  Split  step  methods  (Godunov,  Strang  etc) 

IF  Integrating  factor  methods  (Lawson,  Norsett) 

SL  Sliders  (Fornberg,  Driscoll) 

ETD  Exponential  Time  Differencing  (Cox,  Matthews,  now 
also  Krogstad) 

LG  I  Lie  group  integrators  with  special  actions 
(Munthe-Kaas  and  others). 
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Exact  or  Rational 


Among  all  these  schemes,  some  use  exact  partial  flows, 
others  use  rational  approximants  for  the  unbounded  part. 

Claim 

Real  eigenvalues  of  L  favour  rational  approxmants, 
imaginary  eigenvalues  favour  exact  partial  flows. 


Henceforth  we  consider  only  schemes  which  apply  exact 
partial  flows  (ie.  “exact  up  to  space  discretization”) 
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Assumptions 

In  what  follows,  we  shall  always  assume 

•  Whenever  L  =  0  the  scheme  reduces  to  a  classical 
RK  scheme  for  the  problem  ut  =  N(u ) 

•  Whenever  N(u)  =  0,  the  exact  solution  of  ut  =  Lu 
is  recovered. 


The  classical  Runge-Kutta  scheme  obtained  when  L  =  0 
is  denoted  “The  underlying  RK-scheme” 

Our  favourite  choice  for  underlying  RK  scheme  is  the 
classical  RK4. 
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Classical  Runge-Kutta  4  (RK4C) 

Problem  of  form  ut  =  N(u).  Step  from  t0  to  t0  +  h 


Nj.  = 
N2  = 
N3  = 
N4  = 

ui  = 


N(u0) 

N(u0  +  §fcN4) 
N(u0  +  §ftN2) 
N(u0  +  hN3) 


uo  +  §(Ni  +  2N2  +  2N3  +  N4) 


-2.5 


-1.5 


-0.5  0 
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An  Integrating  Factor  Scheme  (LAW4) 

Lawson  (1967)  derived  the  schemes  by  setting 

v(t )  =  exp  (-tL)u(t) 

which  leads  to  vt  =  N(v)  where  N  =  etL  o  N  o  e~tL. 
Solve  resulting  equation  by  RK4C. 

Ni  =  N(uq) 

N2  =  N(e$L(u0  +  ±hN1)) 

N3  =  N(e*Luo  +  |/iN2) 

N4  =  N(ehLu0  +  e$LhN3) 

U!  =  ehLu0  +  |(ehiN!  +  2e^(N2  +  N3)  +  N4) 
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Lie  Group  Methods  and  the  Affine  Action 


The  scheme  is  based  on  the  affine  Lie  group  action. 
Discrete  case:  Let  G  be  a  matrix  group  over  C  with  Lie 

algebra  g.  Pairs  (M,  b)  e  G  x  CN  act  on  points  in  CN 

( M ,  b)  •  x  =  Mx  +  b 

The  Lie  algebra  consists  of  pairs  (A,  b)  e  g  ix  CN . 
Exponential  map 

—  1 

Exp (t(A,  b))  =  (etA, — —  b). 

Commutator 

[(Ai,  bi),  (A2,  ^2)]  =  [A1A2  —  A2A1,  A162  —  ^4.2^i] • 
Here,  set  g  =  span{X}. 
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An  RK-Munthe-Kaas  Scheme  (RKMK4) 

From  Munthe-Kaas  &  Owren  (1999)  we  derive 
Nx  =  N(u0) 

hL  ,T 

N2  =  N(e  2  Wo  +  1 0o(^)N1) 

Ci  =  hL( N2  -  Ni) 

N3  =  N{ef  uo  +  <M^)(|N2  -  |Ci)) 

N4  =  N(ehLu0  +  MhL)hN3) 

C2  =  h  L{ Ni  -  2N2  +  N4) 

wi  =  ehiw0  +  |</>o(^)(Ni  +  2N2  +  2N3  +  N4  -  Ci  -  |C2) 
where  4>o(z)  = 
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A  Commutator-Free  Lie  Group  Integrator,  Cf4 

Celledoni  et  al.  (2002) 

Ni  =  N(u0) 

hL 

u2  =  e  2  u0  +  f  <M^)Ni 
N2  =  AT(U2) 

hL 

N3  =  N(e  2  u0  +  |0o(^)N2) 

N4  =  N(e^  U2  +  ^0(^)(N3  -  |Ni)) 

hL 

Us  =  e  2  Mo  +  ^L0o(^)(3Ni  +  2N2  +  2N3  -  N4) 

hL  U  U  A 

Ul  =  e  2  Us  +  ^0o(^)(-Ni  +  2N2  +  2N3  +  3N4) 


Applications  of  Lie  group  integrators  and  exponential  schemes  -  p.  23/60 


A  Cox  and  Matthews  Scheme  (C-M4) 

This  scheme  has  the  same  Ni, . . . ,  N4  as  Cf4. 


«i  =  ehAu0+h(f2(hL)N1+2f3(hL)(N2+N3)+f4(hL)N4) 


where 


f2(z)  - 

fs(z)  = 


U(z) 


—  I  —  z  +  ez(4  —  3z  +  z2) 

z 3 

2  +  z  +  ez(-2  +  z) 


—4  —  3z  —  z2  +  ez(4  —  z) 


Derivation  technique:  Unknown! 
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A  unified  format 


By  carefully  studying  all  these  schemes,  one  finds  that 
they  all  fit  into  the  framework 

S 

Nr  =  N(eCrhL  u0  +  hJ2  Q>r(hL)  Nj),  r  =  1, . . . ,  s 

.7=1 

S 

Ul  =ehLu0  +  h^2  br(hL)  Nr. 

r  —  1 


4(z)  =  v  <4’m  zm,  v(z)  =  v 


m  nr 
Z 


rn 


rn 


(a£°),  (/3r,°)  underlying  RK  scheme. 


Applications  of  Lie  group  integrators  and  exponential  schemes  -  p.  25/60 


Order  theory 


Order  conditions  and  B-series  can  be  derived  by 
standard  tools  (rooted  trees). 

T  :  The  set  of  bicolored  rooted  trees  where  each  white 
node  has  at  most  1  child. 

T'  :  Subset  of  T  where  each  white  node  has  precisely 
one  child  (no  white  leaves) 


•  B-series  indexed  by  T. 

•  Order  conditions:  T'  suffices. 


Applications  of  Lie  group  integrators  and  exponential  schemes  -  p.  26/60 


Order  conditions  for  exponential  integrators 

An  exponential  integrator  has  order  p  if 


ui(t)  =  — — ,  for  all  r  g  T'  such  that  |r|  <  p, 


7(r) 


where 


ui(0)  =  Ur(0)  =  1,  1  <  r  <  s, 
u  =  ^/3r’mUr(n)...U  r(rM) 

r 

U r(wqrB+(n,  ■  ■  • ,  T„))  =  £  a^,mUj(ri)  ■  ■  ■  U3-(rM) 

3 
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Number  of  conditions 


Generating  function  for  #  trees  with  q  nodes  in  T' 


M(x)  = 


X 


1  —  X 


,  ,  M(x^)  , 

exp  M[x)  H - 1 - 1- 


The  number  of  order  conditions  for  each  order  1  to  9  is 
1,  2,  5, 13,  37, 108,  332, 1042,  3360. 
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Coefficient  function  spaces 


aJr(z),  br(z )  belong  to  some  function  spaces  we  denote 
Va,  Vb  of  finite  dimension. 

Often,  Va,  Vb  they  are  related  to  the  functions 


Mz)  = 

[  e^~e>>z0kd0, 

Jo 

e  - 

e.g.  Mz)  = 

Z 

Scheme 

Va 

Vb 

Cf4 

Ml),  zMD2 

Ml),  Mz) 

C-M4 

As  Cf4 

Mz ),  Mz),  Mz) 

RKMK4: 

"0- 

O 

fcO|M 

o 

fcO|M 

4>o(z),  zMz) 

Law4: 

1,  ez/2 

1,  ez/2,  ez 
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Assumption  and  bounds 


Let  V  of  dim  K  be  a  function  space  as  above. 
Assumption.  The  map 

/€Vh  (/( 0),  /'( 0), . . . ,  0))  6  R* 

is  injective 

Theorem.  For  any  pth  order  exponential  integrator,  one 
has 


Ka  =  dim  Va  > 


Kb  =  dim  Vb  > 


p  +  1 
2 


Moreover,  the  lower  bound  for  Vb  is  always  attainable  with 
basis  *  •  *  *  *pKb— i 
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Remarks 


•  We  have  no  general  proof  that  lower  bound  for  Ka  is 
sharp.  However,  with  p  =  5  one  can  use  Kb  =  2  with 

•  A  procedure  for  constructing  exponential  integrators 
has  been  developed.  One  starts  with  an  arbitrary 
underlying  scheme  as  well  as  Va,  Vb. 

•  The  really  interesting  part  is  still  ahead:  Choose 
spaces  Va,  Vb  to  deal  with  unbounded  L.  In  the  time 
to  come,  we  focus  in  particular  on  the  NLS. 


Applications  of  Lie  group  integrators  and  exponential  schemes  -  p.  31/60 


Natural  Continuous  Extensions  Zennaro  1986 


Let  (a£°)  (/3r,°)  define  an  underlying  Runge-Kutta 
scheme  of  order  p 

Suppose  that  polynomials  wi(0),...,  ws(0)  of  degree  d 
can  be  found  such  that 

N(to  +  Oh)  :=  J2wr(e)N r 


satisfies 


1 1 


max 


I  N(u(t)) 


G(t)(N(u(t))  - 


-  N(t)\  =  Oih^1) 
N(t))dt  =  G(hp+1 ) 


J  to 

NB!  Requires  aJr(z)  to  be  given. 
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NCEs  continued 


Replace  ut  =  L  u  +  N{u)  by  vt  =  Lv  +  N(t),  solve 
exactly,  and  set  u\  :=  v{h).  Yields  exponential  integrator 
of  order  p  with 

br (z)  =  f  exp((l  —  0)z)w'r (0)  dO 

Jo 

=  (3r,°  +  z  (  exp((l  —  0)z)wr(0 )  d0 

Jo 

In  particular,  these  are  expressed  in  terms  of 

c pk(z )  =  [  exp(z(l  —  0))  0k  d0 ,  k  —  0, 1, . . . 

We  have  rediscovered  the  Cox&  Matthews  schemes. 
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Why  it  works 


(u  —  v)t  =  L(u  —  v)  +  (N(u)  —  N(t )) 

(u  -  v)(ti)  =  e^-^L (N(u(t)  -  N(t ))  dt 

=  0(hp+1) 

thanks  to  the  definition  of  NCEs 
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Theorem  of  Zennaro 


The  following  result  gives  us  a  sharp  lower  bound  for  the 
number  of  </>*.  functions  which  must  be  included 

Any  NCE  satisfies 


<  d  <  min{p,  s*}. 


Moreover,  and  NCE  of  degree  q  always  exists. 


Conclusion.  An  underlying  RK  scheme  of  order  p  can 
always  be  extended  to  an  exponential  integrator  of  order  p 
where 

br(z)  e  span{</>0, . . . ,  Vr 
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The  nonlinear  Schrodinger  equation 

Generally 


iut  =  —A  u  +  (V(cc)  +  A 


u 


2<t 


)  u, 


(t,  x)  eRxRd 


Here 


0  <  cr, 


and  <7  < - 

d- 2 


if  d  >  3. 


•  Potential:  V(cc)  e  Lx  +  L0 c. 

•  1C:  0)  =  «o  6  S  C  J3-1. 

•  Here,  let  d  =  1  and  (#,  t)  6  S'1  x  R. 

•  Usually,  take  a  =  1  (cubic  case). 
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Spectral  Discretisation  in  Space 


Use  2 n  modes  and  set 


71—1 

cfc(t)  =  V  U(^,t)e~imh 


m—  —  n 


leading  to  the  NLS  spectrally  discretised  system 


dc 

d£ 


=  Dc  +  fnoNo5))1(c) 


D  =  diag(-ife2)”"i 


n 


iN(U)e  =  (V(x£)  +  \\U(x£)\2)U(xi). 
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H«P(-,0.1)-'Ph(-,0.1) 


NumExp  1 


Timestep  h 
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A  Simplified  Case 


Let  us 

•  Focus  on  one  scheme,  say  the  Cf4  scheme. 

•  For  analysis,  set  V(x)  =  v  and  A  =  0. 

In  this  case,  the  SDNLS  decouples  into  scalar  equations 

ck  -  akck  +  (3k 

where  ctk  =  —ik2,  and  (3  =  —iv. 

Setting  ak  =  —,  mk  =  e_2hfe2  the  Cf4  scheme  is 

5 

Cl  =  p(mk,ak)c0,  p(m,a )  =  ^  rj(a)mj,  rj  G  n4[o] 

j=o 
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Global  Error  Cf4 


Need  to  find  global  error  at  t  =  T. 

Must  estimate  \p%  —  e%\  where  n  =  T/h 

Pk  =  p(mk ,  ak),  ek  —  exp (-ih(k2  +  v)) 

A  rigorous  analysis  shows  that  up  to  leading  order 


whenever  hk2  <C  1  whereas  for  hk2  1  (and  Ivl  <  \k2) 
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The  Global  Error  for  Decoupled  Case 


h 
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Summing  It  Up 

The  ^2-norm  of  the  global  error  is  found  by  summing  up 


We  may  assume  that  |cg|  <  ^  (holds  in  particular  if  u0  is 

C^S1,  C)).  Assuming  that  N2h  >  1  we  estimate  by 
Euler-Maclaurin’s  formula 


2p— 1 

C  h  4  ,  p  <  8. 
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||«P(-, 0.1) -¥„(•, 0.1) 


Figure  1 


Timestep  h 
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im-,o.i)-'Fh(-,o.i) 


Figure  2 


Global  error,  N=1024,  Initial  condition:  exp(2sin(x)),  Potentials /4||x|p 
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||«r(-,o.i)-'rh(-,o.i) 


Figure  3 


Timestep  h 
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H«P(-,0.1)-'Ph(-,0.1) 


Figure  4 


Timestep  h 
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H«P(-,0.1)-'Ph(-,0.1) 


Figure  5 
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im-,o.i)-'Fh(-,o.i) 


Figure  6 


Global  error,  N=1024,  Initial  condition:  C3-function,  Potential=1/4||x||2 
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H«P(-,0.1)-'Ph(-,0.1) 


Figure  7 
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im-,o.i)-'Fh(-,o.i) 


Figure  8 


Global  error,  N=1024,  Initial  condition:  (^-function,  Potential=1/4||x||2 
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H«P(-,0.1)-'Ph(-,0.1) 


Figure  9 
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im-,o.i)-'Fh(-,o.i) 


Figure  10 


10 


10 
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10 
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10 
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Global  error,  N=1024,  Initial  condition:  (^-function,  Potential=1/4||x||2 
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Our  Work  Plan 


•  Do  this  analysis  for  the  other  schemes 

•  Extend  to  more  general  V{x) 

•  Can  the  analysis  include  the  nonlinear  term? 
Numerical  tests  show  that  it  has  little  impact. 

•  We  want  to  do  the  semiclassical  case 

•  More  space  dimensions,  more  general  boundary 
conditions 

•  Compare  results  with  those  from  geometric 
integrators 
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Long  Time  Integration 


For  systems  with  a  symplectic  structure,  experience 
shows  that  retention  of  such  a  structure  in  the  numerical 
scheme  enhances  the  long  time  quality  of  integration. 
Consider  (  slas  et  al)  the  following  1 D  nonlinear 
Schrodinger  equation 


on  the  circle.  It  is  a  completely  integrable  Hamiltonian 
system  in  the  variables  ( q ,  q *)  with 


q 


4 


L 
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Possible  Strategies 


1 .  Find  a  multisymplectic  structure  Reich  (2000)  and 
discretise  simultaneously  in  space  and  time  to 
preserve  this  structure. 

2.  Introduce  an  integrate  space  discretization 
(Ablowitz/Ladik),  the  resulting  ODE  system  is  not 
canonical.  So  standard  symplectic  integrators  cannot 
be  used,  one  may  then 

(a)  Use  a  Lie-Poisson  type  integrator  based  on  the 
generating  function  technique 

(b)  Apply  (locally)  the  Darboux  transformation  to 
obtain  a  standard  canonical  Hamiltonian  system  to 
which  a  symplectic  scheme  can  be  applied 
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A  Multisymplectic  Structure 


The  above  Schrodinger  equation  can  be  written  in  the 
form 

M zt  +  K zx  =  Vz5(z),  z  G  M4. 

Here  M  and  K  are  skew-symmetric  4  x  4-matrices,  and 

S  :  M4  ->  R. 


The  two  matrices  define  symplectic  structures  w,  n  on 

Rrank(M)  apd  Rrank(K)  _ 

uj(U,  V )  -  VtM U,  k(C7,  V)  -  VtK U. 

In  NLS  we  let  a  =  Re  u,  b  =  Im«,  z  =  (a,  b,  ax,  bx)T  SO 

S{z)  =  \{al  +  bl  +  {a‘>  +  b‘>f). 
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Multisymplectic  Cont’d 

Let  U,  V  be  two  solutions  of  the  variational  equation 

Md  zt  +  Kd  zx  =  ~DzzS{z)dz 
It  easily  follows  that  this  pair  of  solutions  satisfies 


dtu>(U,V)  +  dxn(U,V)  =  0 

the  symplectic  conservation  law. 

The  simplest  multisymplectic  scheme  is  obtained  as  the 
concatenated  midpoint  rule: 
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*0|K 


The  Ablowitz-Ladik  Truncation 


For  simplicity,  set  q  =  u  and  p  =  u*.  The  above  NLS  now 
has  the  formulation 

iqt  =  -qxx  -  %q2p 

-ipt  =  -Pxx  -  Zp2q 


Semidiscretised  version 


iqn  =  - 

-ipn  =  - 


qn—i  qn + 1  2  qn 

( Ax )2 

Pn— 1  “F  Pn-\- 1  2 pn 

(Ax)2 


Pnqn(qn— 1  “F  9n+ 1) 
Pnqn(Pn—  1  “F  Pn+l) 
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AL  Truncation  Cont’d 

It  has  a  noncanonical  Hamiltonian  form 

i  =  P(z)VH(z),  z  -  ( p,q ) 


0 

i 

(Acc)3 


£(*  Pn  ( Qn  —  1  “I”  Qn-\- 1) 

n 


2  ln(l  +  (A x)2qnpn) 


P(z) 


1  +  (A  X^qnPn 

Ax 


N 


n— 1 
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Conclusions — Structure  Preserving  Integrators 


•  The  literature  (Islas  et  al.)  report  excellent  long  time 
behaviour  of  the  presented  approaches. 

•  Which  approach  is  best  depends  on  the  parameters  of 
the  problems  and  the  properties  of  interest. 

•  It  is  admitted  that  these  geometric  integrators  are 
more  expensive  than  classical  ones. 
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